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We present an accurate, first-principles study of the electronic structure and absorption spectrum 
of bulk copper within Density Functional Theory in the Local Density Approximation (DFT-LDA) , 
including the study of intraband transitions. We construct norm-conserving pseudopotentials (PPs) 
including the 3d shell (and optionally the underlying 3s and 3p shells) in the valence, and requiring 
a relatively small plane- waves basis (60 and 140 Rydbergs cutoff, respectively). As a consequence, 
these PPs are strongly non-local, yielding to macroscopically wrong results in the absorption spec- 
trum when momentum matrix elements are computed naively. Our results are compared with exper- 
imental photoemission, absorption and electron energy loss data, and suggest non trivial self-energy 
effects in the quasiparticle spectrum of Cu. 
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I. INTRODUCTION 



Copper has played for a long time a central role in 
the elucidation of the electronic structure of solids. It 
is a relatively inert material, very easy to handle experi- 
mentally; its energy dispersion relations (including spin- 
orbit interaction for some bands) have been measured 
with considerable precision; lifetimes of the band states 
as a function of the distance from Fermi level have been 
determined; surface states have been analyzed in various 
parts of Surface Brillouin zone (for a review see^ ) . From 
the theoretical side, the study of noble metals like Cop- 
per using first-principles methods based on plane-waves 
and ab-initio pseudopotentials (PPs) presents some pe- 
culiar complication with respect to the case of simple 
metals or semiconductors. In fact, in addition to met- 
alicity, which implies the use of an accurate sampling of 
the Brillouin zone in order to describe properly the Fermi 
surface, one must also keep into account the contribution 
of d-electrons to the bonding and to the valence band- 
structure. This means that, within the PP scheme, d 
states cannot be frozen into the core part, but must be 
explicitly included into the valence, yielding a large total 
number of valence electrons (11 for bulk copper). Unfor- 
tunately, a Cu pseudopotential including 3'"'' shell states 
into the valence part is very steep. Hence, when working 
with a plane-wave basis, the usage of a PP of this kind 
may be computationally prohibitive. 

On the other hand, the use of a pseudopotential with- 
out explicit treatment of 3d-electrons has been shown to 
be unreliable^. Hence, first-principles methods based on 
plane- waves have been used only seldom to treat Cu^. 
However, methods have been devised for the construc- 
tion of softer pseudopotentials'^"^, which make the in- 
clusion of the 3'"'' shell in the valence more affordable. 
The price to be paid is that the construction and use of 
such pseudopotentials, which often display a very strong 



1-nonlocality, is a quite delicate matter. In particular, the 
choice of a reference component^ and the transferability 
checks must be done with care. These difficulties are 
more than compensated by the simplicity and elegance 
of the plane-waves formalism in the subsequent calcula- 
tions. In the present work, two possible choices for the Cu 
pseudopotential are explored, i.e: a) including 3s and 3p 
electrons into the frozen core (a quite standard choice), 
and b) including the full 3'"'^ shell in the valence. Fully 
converged Density Functional Theory - Local Density 
Approximation (DFT-LDA) calculations are performed 
at 60 and 140 Ry cutoff, respectively for case a) and b). 

The paper is organized as follows: in Section II we 
give the details of the construction of the pseudopoten- 
tials used throughout this work ; in Sections HI and IV 
we present the ground state properties and bandstruc- 
ture, respectively, obtained with the different pseudopo- 
tentials; finally, in Section V, we compare the theoretical 
absorption and electron energy loss spectra with the ex- 
perimental data, including the effects of Local Fields and 
intraband transitions. 



II. PSEUDOPOTENTIAL GENERATION 

The Cu atom has the ground-state electronic config- 
uration [kv]M^°As^ (or [Ne]3s23p63di°4si). 4s and 3d 
eigenvalues are separated, in DFT-LDA, by less than 0.5 
eV. It is then quite obvious that freezing all states but the 
4s one into the atomic core (i.e., to neglect the polariza- 
tion of the 3d electrons) cannot yield a good, transferable 
pseudopotential. On the other hand, inclusion of the 3d 
electrons into the valence (we call this a "3d" pseudopo- 
tential) yields a much slower convergence of plane-wave 
expansions, due to the steepness of the d-component of 
the PP. Unfortunately, the spatial superposition between 
the 3d and 3s or 3p states is quite large, despite the large 
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(~ 70 eV) energy separation. Hence, an even more con- 
servative and secure choice for the PP is to include all 3s, 
3p and 3d electrons into the valence: in fact, 2p and 3s 
states are well separated, both spatially and energetically 
(~ 800 eV). This latter choice gives rise to a PP which 
is even harder than the "3d" one, and which will be re- 
ferred to as a "3s" pseudopotential, yielding 19 valence 
electrons per atom. 

Using the traditional pseudopotential generation 
scheme proposed by Bachelet, Hamann and Schliiter 
(BHS)^ yields PP whose d component converges very 
slowly in Fourier space, requiring to work at an energy 
cutoff of 200 rydbergs or more. However, PPs which con- 
verge at less than 100 Rydbergs can be constructed by 
using specially devised schemes as those of refs.'^"^. We 
choose to restrict to the class of norm conserving pseu- 
dopotentials, in order to avoid the additional numerical 
complications arising from the charge-state dependence 
of the PP of ref.^. Moreover, our choices in generating the 
PP are dictated not only by the need of a fast convergence 
of the PP in Fourier space, but also by that of optimizing 
the PP accuracy and transferability, a non-trivial task 
when also 3s and 3p states are included into the valence. 
In particular, the Hamann scheme* (which does not try 
to optimize the Fourier space convergence at all) turned 
out to yield much more accurate and transferrable norm- 
conserving PPs, particularly when the whole 3'"'* shell is 
included into the valence. 

Hence, we strictly follow the Hamann procedure (de- 
scribed in the appendix of ref.*) whenever it is possible, 
i.e. in all cases except for the d components. For the lat- 
ter, we follow the TrouUier-Martins^ scheme, which al- 
lows us to reduce significantly the number of plane-waves 
requested for convergence without loosing too much in 
transferability. In all other cases (i.e. for the s and p 
components of both "3s" and "3d" pseudopotentials) we 
found the Hamann procedure to be more convenient, even 
at the cost of a slower Fourier space convergence, since 
the TM one yielded significantly worse results, and/or 
ghost states^ when the PP were used in the Kleinman- 
Bylander form^°. 

Another delicate point is the choice of a reference com- 
ponent, i.e. of a PP angular momentum component 
which is taken to be valid for every 1> 3. Often, the 
1=2 component is chosen as a reference, simply because 
this makes calculations easier. This choice is sometimes 
lacking a physical justification, and can be dangerous, as 
it has been shown by some of us in the case of Sn^^ . In the 
present case the 1=2 reference had to be avoided anyway, 
since it yielded a much worse transferability than the 1=0 
or 1=1 choices, and sometimes gave rise to ghost states 
in the KB form. 

Several trials and tests with different cutoff radii have 
been done, in order to optimize the PP transferabil- 
ity without increasing too much the number of plane 
waves requested for convergence. Transferability tests in- 
cluded both the plot of logarithmic derivatives, and the 
explicit calculation of pseudoatom eigenvalues in some 



excited electronic configuration (both neutral and posi- 
tively charged). 

Finally, since 3s and 3p states in the solid preserve their 
atomic configuration better than 3d ones, their explicit 
inclusion into the valence can sometimes be avoided, by 
considering only the effects of the non-linearity of the 
exchange-correlation potential^^. Hence, a third pseu- 
dopotential, with frozen 3s and 3p electrons but includ- 
ing non-linear core-corrections has also been considered, 
and will be referred to as the "3d-FNLCC" PP. In the 
latter, the core charge is represented by the true one for 
r> 0.5 Bohrs, and by a Gaussian model charge for r< 0.5 
Bohrs. 

The resulting optimal cutoff radii and reference 1- 

components, chosen for our "3d" , "3s" , and "3d-|-NLCC" 
pseudopotentials, are given in Table I. 
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FIG. 1.: Calculated total energy vs. lattice constant for bulk 
Cu. Panel (1): effect of the kinetic energy cutoff at fixed number 
of k-points (28). Panel (2): effect of the IBZ sampling at fixed 
Ecut (75 Ry.); Panel (3): effect of nonlinear core-corrections (sec 
text). The reported values for the equilibrium lattice constant and 
bulk modulus have been obtained from Murnaghan fits (continuous 
curves). The experimental values are a=6.822 Bohrs and B=1.827 
MBar^s. 



III. GROUND STATE PROPERTIES 

Our first step is a self-consistent ground-state calcu- 
lation, performed by minimizing the DFT-LDA energy 
functional with a Car-Parrinello method"'^^, in a stan- 
dard plane-wave basis. The Ceperley-Alder^^ exchange- 
correlation energy and potential, as parametrized by 
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Perdew and Zunger-"^^, have been used (test calculations 
with the Hedin-Lundqvist form^^ have also been per- 
formed: see below). All our pseudopotentials are used 
within the fully-separable Kleinman-Bylander scheme^°, 
after checking that no ghost-states was present^. The ir- 
reducible wedge of the Brillouin zone (IBZ) was sampled 
with the use of Monkhorst and Pack (MP) sets^^ of 
k-points. A very small fictitious electronic temperature 
(equal to ~ 10 Kelvin) was used in order to accelerate 
the convergence of the calculated Fermi surface. Conver- 
gence with respect to both the k-points sampling and the 
kinetic energy cutoff has been checked extensively: Fig. 1 
and Table II show the results obtained for the "3d" pseu- 
dopotential. Ec„t = 60 Ry. and — 28 appear to give 
well-converged, satisfactory results, with lattice constant 
ao and Bulk modulus Bq within 1.4% of the experimental 
values. In the calculation with NLCC, the energy cutoff 
had to be increased to 100 Ry, in order to describe prop- 
erly the core charge. NLCC reduce the underestimation 
of the experimental lattice constant to 0.5 %, but they 
are found to be almost ininfluent on the bandstructure, 
as well as on the resulting spectra (see below). Finally, 
calculations with the "3s" pseudopotential, much deeper 
than the "3d" one, required an energy cutoff of 140 Ry. 
Also in this case, the effect on the LDA bandstructure 
and spectra in the energetic region of interest are found 
to be very small (see below). 

In summary, as long as only the LDA 4s and 3d band- 
structure is concerned, the effect of the inclusion of 3s 
and 3p in the valence, (as well as the use of NLCC), are 
mainly confined to a change of the equilibrium lattice 
constant, which induces an indirect effect on the band- 
structure energies^*. 

bS' = 60 Ry 
AE = 0.014 eV 

1 

AB = 0.15 eV ,,HL 



without NLCC 

1 

A£ = 0.214 eV 

with NLCC 



PP without 3s3p 

! 

A£ = 0.2 eV 

PP with 3s3p 



FIG. 2 . : Summary of the effects of different computational de- 
tails on the LDA bandstructure. The energy differences reported 
are the maximum ones, and correspond to the bottom valence band 
at the r point (Fi). They decrease gradually to zero at the Fermi 
level. 



IV. LDA BAND STRUCTURE 

Fig. 2 summarizes the maximum relative energy differ- 
ences induced by different computational details on the 
LDA bandstructure; the value is taken at the bottom va- 
lence in the F point (Fi) and decreases gradually to zero 
at the Fermi level. 

By changing E^ut from 60 Ry to 75 Ry the 
band-structure remains practically identical (changes 
are less then 0.02 meV). At the 60 Ry cut- 
off, we compare the band structure calculated with 
the Ceperley-Alder parametrization of the exchange- 
correlation potential^'*' with that obtained using the 
Hedin-Lunqvist parametrization^^; in the latter case we 
find a maximum energy shift of about 0.15 eV. 

The inclusion of NLCC yields, instead, a maximum 
upward shift of about 0.21 eV with respect to the case 
without NLCC, mainly due to the change in the equilib- 
rium lattice constant. To study more deeply the effects 
of core polarization, we also performed a band-structure 
calculation with the "3s" pseudopotential (see Section 
II), where the 3s and 3p core level relaxation is fully in- 
cluded in the selfconsistent run. Also in this case the 
maximum bandshift with respect to the "3d" pseudopo- 
tential is limited to about 0.2 eV. Hence, our results ob- 
tained with the "3d" pseudopotential at i?ctit=60Ry. and 
using 28 MP k-points in the determination of the self- 
consistent charge density can be considered to represent 
the converged LDA bandstructure of bulk copper. Theo- 
retical results are compared with the experimental pho- 
toemission data in Fig. 3. 




L r X K r 

FIG. 3.: Bulk copper DFT-LDA bandstructure (full-line), com- 
pared with photoemission Data (points) from Rcf. ^. 

At difference with the case of semiconductors, the dis- 
agreement between theory and experiment is far from 
being limited to a rigid shift of the Kohn-Sham occupied 
eigenvalues with respect to the empty ones. In fact, as 
also summarized in Table III, the widths of the d bands 
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are systematically overestimated, a well-known failure of 
LDA when applied to transition and noble metals^^. 



V. ABSORPTION SPECTRUM 

The absorption spectrum is given by the imaginary 
part of the macroscopic dielectric function: 



CM (w) 



EO ('^) ' 



where 



4n 



('^) = 1 + 1™ ■r-F^XG=o G'=o (q, 



q^o |q 



(1) 



(2) 



Xg g' (q, is the reducible polarization, solution of the 
equation 



Xg G' (q, w) = Xg G' (q> + 

47r 



Yl G" (q^ 



77T2XG"G'(q,w). (3) 



G' 



G"|2 



Neglecting (for the moment) intraband transitions, 
Xg g" (q^ is given by: 

XGG'(q,'^) = 2 / 77^3 

^ Jbz (ztt) 

J2 (^'k - q|e-*('i+°>'"|^k)(nk|e*(^+^')-'|n'k - q) 

e?L(n,n',k,q,w), (4) 

with 

Gl (n, n', k, q. CO) = fn' (k - q) (2 - /„ (k)) 

1 



w + e„' (k - q) - e„ (k) + irj 
1 



e„ (k) - £„' (k - q) - 277 



(5) 



with < /„ (k) < 2 representing the occupation num- 
ber summed over spin components. The sums over k are 
transformed to integrals over the BZ, and the latter arc 
evaluated by summing over large sets of random points 
contained in the whole BZ. Fully converged calculations 
with a small broadening require a very large number of 
k-points^°; in the present work, the broadening used (and 
the corresponding number of k-points which was found 
to be sufficient to ensure convergence) is specified explic- 
itly for each one of the reported spectra. 

The simplest approach to the calculation of the ab- 
sorption spectrum neglects the full inversion of Eq. (3) 
(i.e., neglects Local Field Effects), and assumes 



Cm (w) 1 - lim |— |2XG=o G'=o (q, w) 
q^o |q|^ 



(6) 



The q ^ limit for the oscillator strengths appearing in 
Eq. (4) is calculated in the transversal gauge^^, within 
first order perturbation theory^^ 

lim(n'k-q|e-*'in"k) = 

where (/)„k (r) are the Bloch functions. Due to the non 
local character of the norm-conserving pseudopotentials, 
the well-known relation between [r, H] and the momen- 
tum operator 



[r,H]=p. 
must be substituted by: 

[r,H]=p+[r,VNL]. 



(8) 



(9) 



The second term of the rhs of Eq. (9), which in simple 
metals and in many semiconductors is small (and often 
neglected in practical calculations), becomes extremely 
important in the case of copper due to the large nonlo- 
cality of the PP. 




FIG. 4.: Full line: imaginary part of the macroscopic dielectric 
function of Cu without Local-Field effects and without including 
the non-local pseudopotential comiimtator [Eq. (9)], compared 
with experimental data (squares) from'*". Dot, dashed and long- 
dashed lines correspond to the functions defined in Eq. (10) with 
n' = 3, n' = 4 and n' = 5 respectively. All theoretical curves are 
computed with Nj, = 15,386, and a Gaussian broadening of 0.15 
eV (see text). 



This is demonstrated in Fig. 4 where the imaginary 

part of sm is calculated assuming the validity of Eq. 
(8) The experimental absorption spectrum is severely un- 
derestimated between the offset of interband transitions 
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(« 1.74 eV) and 5 eV. A better analysis of this behavior 
can be performed by plotting the quantity 



EkeL(//,/*',k.O,u;) 



(10) 



with n = 6 and n' = 3,4,5 (see Fig. 4). This quantity 
is the joint density of states (JDOS) (divided by w^) for 
transitions between the sixth band (d-like) and the third 
to fifth bands (sp-like), and corresponds to assuming an 
oscillator strength equal to one in Eq. (4). 




<o [eV] 



FIG. 5.: Full line: imaginary part of the macroscopic dielectric 
function of Cu without Local-Field effects and including the non- 
local pseudopotential commutator [Eq. (9)], compared with exper- 
imental data (squares)^". Nj. and broadening as in Fig. 4. 



We see that transitions in the energy range of interest 
(1.8 to 5.0 cV) exist, but they arc strongly suppressed due 
to the small values of the corresponding matrix elements 
of p. Using both terms of Eq. (9), instead, we obtain the 
spectrum shown in Fig. 5, comparing much better with 
experiments. The large influence of the non-local pseu- 
dopotential commutator on d ^ s/p optical transitions 
can be understood by writing explicity the contribution 
of the second term of Eq. (9): 



{(t>{s/p)k\ [r, Vnl] \(t>id)k} = jj drdr'cl)l^/p^j^ (r) 

[r\4^ (r, r') - (r, r') r'] (r') , (11) 



where (r, r') is the 1-orbital component of the pseu- 
dopotential, and (/'(s/p)k (r)) 'P(d)k (i") ^tre the s/p like and 
d like Bloch functions. Now we can approximate the sum 



in Eq. (11) with the leading terms, to obtain: 

{4>is/p)k\[Tc,VNL]\(l){d)k} ^ jj drdr' 

}(sMi.i^)^VNL (r,r').^(d)k(r')- 
'^(Vp)k (r) (V^L (r, r') + V^^ (r, r')) r>(,)k (r'] 



(12) 



In the case of copper the d components of the pseudopo- 
tential differ from the p component near the origin by 
about 20 Hartrees, and this explain the strong influence 
on (w). 

Despite the strong improvement of the agreement with 
experiment obtained in Fig. 5, the theoretical curve ex- 
hibits an amplitude overstimation of about 20% with re- 
spect to the experimental absorption spectrum. 
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FIG. 6.: Effect of the use of a pseudopotential including the 3s 
and 3p shells in the valence on the calculated e (ui) for bulk Cu 
(full Une : "3d" PP; dotted Une: "3s" PP. See text). Due to the 
large number of plane-waves required by the "3s" PP, the compar- 
ison is done using a small number of k-points and a relatively large 
Gaussian broadening (N^ = 3,000, 7 = 0.2 eV). Local-Field effects 
axe neglected, and Eq. (9) is used. 

This drawback must be analyzed taking into account 
both the physical approximations involved in our theo- 
retical approach, and the possible residual errors due to 
the PP scheme. Concerning the first ones, the most im- 
portant point is the neglection of self-energy effects in 
the bandstructure calculation, and of cxcitonic effects 
in the absorption process. Concerning the PP scheme, 
instead, a possible reason for the overstimation of the 
spectrum intensity could be related to the use of pscudo 
wavefunctions instead of the all-electrons ones in Eq. (4). 
This effect has been studied in atoms by the authors of 
Ref.^^: they found that PP calculations, even when the 
second term of rhs of Eq. (9) is correctly taken into ac- 
count, can be affected by a small residual error due to the 
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difference between all-electron wavefunctions and pseudo 
wavefunctions inside the core region. In the case of the 
Cu atoms, 3d Ap transitions were found to yield a 
matrix element which was too large by about 10%^^. 

The effect of this overstimation of the — !■ 4p intra- 
atomic optical matrix elements on the calculated bulk 
spectrum is, however, not obvious. To clarify this point, 
we have performed an accurate comparison of the results 
obtained by using the "3d" and "3s" pseudopotcntials. 
both in case of the bulk crystal and for the isolated Cu 
atom. For the latter case, our results are summarized 
in Table IV: while our "3d" pseudopotential gives about 
the same results as those obtained in Ref.^"^, the "3s" PP 
reduces the error on the intra-atomic optical matrix el- 
ements to less than 1.5 %. However in the case of the 
bulk crystal the amplitude of our calculated spectrum 
does not change appreciably when results obtained with 
"3d" and "3s" pseudopotcntials are compared (see Fig. 
6). This suggests that the overall intensity overestima- 
tion cannot be ascribed to the use of PP wavefunctions 
in Eq. (4). 

The reason for this discrepancy between theory and 

experiment should hence be searched for within the phys- 
ical approximations made, such as the neglection of self- 
energy effects that are currently under investigation and 
will be the subject of a forthcoming paper^'*. 



A. Intraband transitions 

A peculiar characteristic of metals is the intra-band 
contribution to the dielectric function, which, neglecting 
local field effects, is given by: 



^intra (w) = 1 - lini 



477 



q-0 |q|- 



E 



BZ (27r)^ 



with the n-sum restricted only to semi-occupied bands. 

In the case of copper only the sixth band contributes 
to Eq. (13). Indicating this band with np and using the 
relation 

U (k - q) - In, (k) = (k - q) - fn, (k)) 

{e (k - q) - /„, (k)] - e [/„, (k) - (k - q)]}, 

(14) 

time-reversal symmetry allows one to rewrite Eq. (13) as 



f Stt /■ d^k 
\ |qP Jbz (27r)^ 



^intra (w) = 1 - lim 

q^o I |q|^ Jj^y f27r) 
(/„, (k - q) - (k)) e [fn, (k - q) - (k)] 
Knkle^^-lnk - q)|2 (e„, (k) - e„, (k - q)) ' 



- (e„p (k) - e„p (k - q)) 



, (15) 



In the small q limit, Eq. (15) yields the well-known 
Drude contribution to the dielectric function 



e;„,.„(a;)^ 1-^ + (,2) 



(16) 



with 



ujj-, = lim 
q 



im /— / '^^^ 
™ \ |qP Jbz (27r)^ 



(/„, (k - q) - (k)) e [/„, (k - q) - (k)] 
|(nk|e*^--|nk- q)|2 (e„, (k) - (k- q))} . (17) 

Since the theta function in Eq. (17) limits strongly the 
region of BZ that contributes to the integral, the k-space 
sampling and the modulus of the chosen q vector used in 
the numerical evaluation of Eq. (17) become two critical 
convergence parameters: for a small |q| very few k-points 
will satisfy the condition [/„^ (k — q) — fnp (k)] ^ 0. 
In practice, |q| must be small enough to reproduce the 
q ^ limit in Eq. (17), but large enough to allow a 
suitable number of k-points to contribute to the sum. 
To obtain a well converged cj|j we found it necessary to 
use « 16, 000 random k-points in a region of the BZ such 
that tnp (k) is contained within epermi ±0.1 eV. The 
value used was |q| = 0.005 a.u.. A fictious electronic 
temperature was introduced to smear out the Fermi sur- 
face, increasing the number of (k, k — q) pairs giving non 
zero contributions to Eq. (17). In Table V we present 
our results as a function of the fictious electronic temper- 
ature. The optical mass, defined as 



rriopt = 



with 



10.8 eV (18) 



with Vt direct lattice cell volume, converges to 1.36, a 
value in good agreement with the experiments^. 



B. Local Held effects 



Local-Field effects are accounted for when the macro- 
scopic dielectric function is computed according to Eq. 
(1-3), i.e. by explicitly obtaining xg G' (qji^) as 



XG G' (q, w) = Xg g' (q> 

47r 



|q + G| 



rXG G' (q^t^) 



, (19) 



We divide explicitly Xq q, (q, w) into intraband and in- 
terband contributions: 

x°G G' (q. '^) = x'g G' (q, ^) + x'g G' (q, , (20) 
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where the interband part is given by Eq. (4), while for 
the intraband contribution wo have to evaluate: 



1 t fl^\< 

^intra , A ^ / 

Xg g' IQ' - 1^ 

^ JBZ ZTT 



BZ (2 

(nj.k - q|e-*('i+^)-""|«Fk)(nFk|e<'i+^')-'|n^k - q) 

("-F,"-F,k,q,tj) . (21) 

with Ql^ {n, n', k, q, w) defined in Eq. (5). For G , G' ^ 
the oscillator strength is calculated using Fast Fourier 
Transforms (FFT) 

(nk - q|e-^('i+°)-'-|n'k) = (u„k|e-^°-nu„'k) + O (q) , 

(22) 

while for the G = G' = element of ^^^^^^ we use Eq. 7. 




FIG. 7.: Real and imaginary part of SMi^) for bulk Cu with (full 

line) and without (dotted line) inelusion of Loeal-Field effeets, com- 
pared with experimental data (squares) from Ref.-"'. Theoretical 
spectra are computed with = 110 (in the irreducible wedge of 
BZ), and a Lorentzian broadening of 0.4 eV (see text). 

The method presented in Section VA for the cal- 
culation of lod could, in principle, be extended to 
Xgg''(^''^)- Unfortunately the explicit calculation of 
^tntra ^j-^^ j^Q^ Q/-J pg^jj-g jg computationally pro- 

hibitive, due to the large number of k-points required 
to reach convergence. To overcome these difficulties, 
we have evaluated Eq. (21) on a limited number of 
k-points (a MP grid of 110 points in the irreducible 
wedge). (rii?, ni?, k, q, w) is different from zero only 
for k very close to the Fermi Surface (which, in the 
grid, coincides with one particular point k^) where 
fnp (k — q) (2 — fnp (k)) 7^ 0. The oscillator strengths 
can be considered to be almost constant in the vicinity 
of the Fermi Surface. The same ansatz cannot be applied 
to Q^^^ {np, np, k, q, cj); however we can use the property 
that near the Fermi surface the metallic band dispersion 
of copper is well approximated by a free-metal one, and 



make the following assumption for Xq (q, w): 

xSG'(q,'^)« 

^(nj.RkF - q|e-^('»+*^)-""|nFRkF) 



1 

n: 



I R 



(nFRkF|e'(''+^')"-'|nirRkF - q)l ttq (q,a;) , (23) 

where R is one of the Ns symmetry operations not in the 
point group of hp- ttq (q, w) is the non-interacting polar- 
ization calculated for a jellium model^^ with a density Uei 
yielding a classic plasma frequency cOp = y/Ann^i = 9.27 
eV, corresponding to the value of cod calculated in Sec- 
tion V A. 



TTO (q, w) 



27r2|q| Jo 



xdx 



CO - k'p |q| +iri uj - k^p |q| 



IT] 



where 



and 



OJ 



kp'' = {rn^nei) 



1/3 



(24) 



(25) 
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FIG. 8.: Imaginary part of the inverse macroscopic dielectric func- 
tion of Cu including Local-Field effects and the Drude contribution 
(full line). Dashed line: results without Local-Field effects; Long- 
short dashes: with Loeal-Field effects, but without the Drude con- 
tribution. All theoretical spectra are computed with the same 
of Fig. 7 but with a Lorentzian broadening of 0.9 eV. Points are 
the EELS data from Ref.^o. 

In Fig. 7 we compare our results for £jv^(w) with and 
without Local Field Effects. The differences are small, 
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consistently with the fact that, as expected, no large 
LFE are present in a metal. However Local Field effects 
are more important on the EELS spectrum, as shown 
in Fig. 8: in the high energy region the full inversion 
of Xg g' (q, w) matrix corrects an overstimation of the 
intensity. The inclusion of intraband transitions, on the 
other hand, turns out to be necessary not only to describe 
correctly the behaviour of e at low frequencies, but also 
to improve the agreement of the calculated e~^{w) with 
the experimental EELS data (Fig.8). 



VI. CONCLUSIONS 

Several conclusions can be drawn from the presented 
results. First, we have shown that fully converged DFT- 
LDA calculations for bulk copper using plane waves can 
be performed without requiring an exceedingly large ba- 
sis set, by using soft norm-conserving pseudopotentials 
including the 3d electrons in the valence. 

The strong non-locality implied by such pseudopoten- 
tials reflects itself in a very large contribution of the com- 
mutator between the PP and the position operator in the 
calculation of the matrix elements entering the dielectric 
function. 

Similar calculations including also the 3s and 3p shells 
in the valence are also feasible using plane waves, but 
show very little difference on the bandstructure and on 
the optical matrix elements not directly involving those 
shells. 

At difference with the case of simple semiconductors, 
the disagreement between the LDA theoretical band- 
structure and the experimental one cannot be corrected 
by a rigid shift of the Kohn-Sham eigenvalues: this sug- 
gests the importance of many body effects in the quasi- 
particle spectrum of Cu. These effects should be taken 
into account through a calculation including self-energy 
contributions^^. By contrast, the computed optical mass, 
which involves only intraband transitions at the Fermi 
energy, is in very good agreement with the experimental 
data. Also the theoretical absorption spectrum, where 
only LDA eigenvalues within a few eV around the Fermi 
energy are important, displays a quite satisfactory agree- 
ment with the experimental one. The overall overestima- 
tion of about 20 percent in the amplitude of the main 
peaks cannot be attributed to the effect of wavefunc- 
tion pseudization inside the core region, and should be 
ascribed to the neglection of self-energy (and excitonic) 
effects. Finally, local field effects on the macroscopic di- 
electric function turn out to be negligible, with or without 
the inclusion of the Drude term. 
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Psi'u<li)potoiitia.l 


r., ^Bulirj 


r^, [Bi)lir] 


r,/ [Biilir] 


RcfiTi'iKx; 1 coiiiptjiiciit 


"3d" 


1.19 (H) 


1.19 (H) 


2.08 (TM) 


s 


"3d+NLCC" 


1.10 (H) 


1.19 (H) 


2.08 (TM) 


s 


"3s" 


0.49 (H) 


0.60 (H) 


1.19 (TM) 


p 



TABLE I. Cu Pseudopotential cutoff radii and reference components. (H) and (TM) stand for Hamann (ref *) and Troul- 
liers-Martins (ref.®) schemes, respectively. All our pseudopotentials are norm-conserving in the sense of Bachelet, Hamann 
and Schliiter'^, and have been produced using a pubhcly available Fortran code^*. They are also available in numerical form 
upon request by e-mail to: onida@roma2.infn.it 



V\\ cutoff (Ry) 


Xk 


flu [B(jlirj 


Bu [MBar] 


75 


10 


6.77 


1.73 


75 


28 


6.73 


1.79 


75 


60 


6.74 


1.83 


60 


28 


6.73 


1.83 


100-l-NLCC 


28 


6.79 


1.71 


Expt. 




6.82 


1.83 


TABLE II. Convergence of the calculated ground-state properties of bulk Cu. 






Experiment^ 


Present work 




ri2 


-2.78 


-2.33 


positions of d-bands 


Xs 


-2.01 


-1.46 






-2.25 


-1.69 




ri2 — r26' 


0.81 


0.91 




X5 — Xz 


2.79 


3.23 


widths of d-bands 


X5 — Xi 


3.17 


3.70 




Li - Li 


1.37 


1.58 




Li - Li 


2.91 


3.72 


positions of sp-bands 


Fi 


-8.60 


-9.85 




L21 


-0.85 


-1.12 


L-gap 


L\ — L2' 


4.95 


4.21 


TABLE III. Comparison of Cu band widths and energy position 


with experimental values^ at high-symmetry points. All 


energies in eV. 








Optical transitions 


AE 


"3d" PP 


"3s" PP 


4,s Ap 


1.7.32 


1.738 


1.713 


:',(/ — ip 


().ll)() 


I). 15:! 


I). 112 



TABLE IV. Optical matrix elements for the Cu atom, calculated in the All-electrons (AE) scheme, and with the pseudopo- 
tentials used in the present work. Values are in atomic units. 
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Tel [eV] Number of intraband transitions ujd [ev] rriopt 

0.1 14,617 9^39 02" 

0.01 2,456 9.41 1.32 

0.001 534 9.25 1.36 

0.0001 310 9.28 1.36 

0.00001 296 9.27 1.36 



TABLE V. Cu Drude plasma frequency and optical mass values obtained using different fictious electronic temperatures; 
the experimental value for niopt is 1.35^^. 
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